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Abstract 

This  paper  describes  the  optimization  of  a 
piezoelectric-driven  synthetic  jet  actuator  based  on  a 
Lumped  Element  Modeling  (LEM).  To  simplify  the 
problem,  this  papers  splits  the  optimization  problem 
into  two  parts.  First,  a  constrained  optimization  of 
the  cavity  volume  and  orifice  dimensions  of  two 
baseline  synthetic  jets,  each  with  a  given 
piezoelectric  diaphragm,  is  conducted  using  two 
different  objective  functions.  One  seeks  to  improve 
the  centerline  output  velocity  over  a  broad  frequency 
range,  and  the  other  maximizes  the  centerline 
velocity  at  a  prescribed  resonant  frequency  of  the 
device.  Significant  improvements  are  achieved  using 
both  objective  functions  for  both  synthetic  jets. 
Second,  the  two  baseline  piezoelectric  diaphragms 
have  been  optimized  using  two  configurations.  One 
uses  the  standard  inner-disc  piezoceramic  patch 
bonded  to  a  metal  shim,  while  the  other  employs  an 
outer  piezoceramic  ring.  In  each  case,  the  objective 
is  to  maximize  the  achievable  volume  displacement 
of  the  diaphragm  at  the  coercive  electric  field 
strength  of  the  piezoceramic,  while  the  natural 
frequency  of  the  piezoelectric  diaphragm  is 
constrained  to  be  greater  than  or  equal  to  the  baseline 
designs.  Both  configurations  yield  modest  (~5%) 
improvements  for  one  diaphragm  and  significant 
improvements  for  the  other  diaphragm  (>50%). 

Nomenclature 

a0  orifice  radius  (mm) 
c0  air  speed  of  sound  (m/s) 

CaC  cavity  acoustic  compliance  =  V0  / p0c0 2 
(s2'm4/kg) 
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CaD  diaphragm  short-circuit  acoustic  compliance 
=  AV/P\k  __0  (s2m4/kg) 

CEF  piezoceramic  electrical  free  capacitance  (F) 
da  effective  acoustic  piezoelectric  coefficient 
=  AV/Vac  |p=0  (m3/V) 

E f  max  maximum  electric  field  (V  /m) 

/  frequency  (Hz) 

f H  Helmholtz  frequency 

=  l/W(M^+M^)Cac  (Hz) 
fD  short-circuit  diaphragm  resonant  frequency 
-l/2^VMoDCaC  (Hz) 

fx ,  f2  synthetic  jet  lowest  and  highest  resonant 
frequencies,  respectively  (Hz) 

Kd  nondimensional  orifice  jet  dump  loss 
coefficient 

L  orifice  length  (mm) 

MaD  diaphragm  acoustic  mass  (kg/m4) 

=  -^r]pA[w(r)J  rdr 

MaN  orifice  acoustic  mass  (kg/m4) 

=  ApQL/?)7ral  for  Poiseuille  flow 
MaRad  orifice  acoustic  radiation  mass  =  %p0/37T2a0 
(kg/m4) 

P  differential  pressure  on  the  diaphragm  (N/m2) 

q  electric  charge  stored  on  the  piezoelectric  (C) 

Qc  volume  flow  rate  through  the  cavity  (m3/s) 

Qout  volume  flow  rate  through  the  orifice  (m3/s) 

Q  volume  flow  rate  displaced  by  the  diaphragm 
=  Qc+Qout  (m3/s) 

RaD  diaphragm  acoustic  resistance 
=  2 C^MJCZ  (kg/sm4) 

RaN  viscous  orifice  acoustic  resistance  (kg/s  m4) 

=  8 F/JqL/^4  for  Poiseuille  flow 
Ra0  nonlinear  orifice  acoustic  resistance 
=  0.5 KDp0Qoul/7T2 a*  (kg/sm4) 

Rx  piezoceramic  radius  (mm) 


R2  shim  radius  (mm) 

s  Laplace  variable  =  jco 

t  piezoceramic  thickness  (mm) 

ts  shim  thickness  (mm) 

u0  centerline  orifice  velocity  (m/s) 

Vac  input  ac  voltage  (V) 

Vac  max  maximum  applied  voltage  =tpx  £/max  (V) 

V0  cavity  volume  (mm3) 

w(r)  transverse  displacement  of  the  diaphragm  (m) 

AV  volume  displaced  by  the  diaphragm 

Rl 

=  J  2 7Trw(r)dr  (m3) 

o 

(j)a  electroacoustic  turns  ratio  of  the  piezoceramic 

diaphragm  =  \da/CaD\  (Pa/ V) 
v  air  kinematic  viscosity  (m2/s) 

pA  area  density  (kg/m2) 

Pq  air  density  (kg/m3) 

co  radian  frequency  =  2nf  (rad/s) 

£  empirical  diaphragm  damping  coefficient 

Introduction 

Synthetic  or  “zero-net  mass  flux”  jets1  are 
commonly  used  as  flow-control  actuators  in  a  wide 
spectrum  of  applications  including  jet  vectoring,2 
separation  control,3’4  and  boundary  layer  control.5’6 
The  performance  specifications  of  any  actuator  are 
quantified  in  terms  of  an  exhaustive  list  of  parameters 
such  as  bandwidth,  control  authority,  etc.  Flow- 
control  applications  require  a  known  actuator 
frequency  response  function  that  relates  the  input 
voltage  to  the  output  property  of  interest  (e.g., 
maximum  velocity,  volumetric  flow  rate,  momentum 
flux,  etc.).  Clearly,  the  required  performance  metrics 
are  application  specific,  and  methods  are  needed  to 
achieve  the  optimal  design  of  these  devices.  Design 
and  optimization  studies  have  been  conducted  for 
piezoelectric  cantilever-type  flow  control  actuators, 
but  the  modeling  issues  are  simpler  compared  to 
synthetic  jets.7  Specifically,  the  cavity  and  orifice 
configuration  of  a  synthetic  jet  significantly 
complicates  the  overall  system  dynamics.  Recently, 
LEM  has  been  combined  with  equivalent  circuit 
representations  to  estimate  the  nonlinear  dynamic 
response  of  a  synthetic  jet  as  a  function  of  device 
dimensions  and  material  properties.8’9  These  models 
have  provided  good  agreement  between  predicted 
and  measured  frequency  response  functions  and  thus 
are  suitable  for  use  as  design  tools. 

The  purpose  of  this  paper  is  to  leverage  LEM  in 
the  optimization  of  piezoelectric-driven  synthetic  jet 


actuators.  The  following  section  briefly  reviews  the 
lumped  element  model  and  corresponding  equivalent 
circuit  of  a  synthetic  jet  from  Gallas  et  al.8  and 
discusses  the  basic  dynamic  behavior.  The  resulting 
system  model  is  then  employed  in  an  optimization 
scheme.  For  the  current  study,  the  optimization 
problem  is  decoupled  into  two  parts.  First,  the 
piezoelectric  diaphragm  is  fixed  while  the  cavity 
volume  and  orifice  dimensions  are  varied.  Second, 
the  piezoelectric  diaphragm  is  optimized  while  the 
other  design  parameters  are  fixed.  Discussions  on  the 
chosen  cost  functions  and  accompanying  constraints 
employed  in  the  optimization  and  the  corresponding 
results  are  presented.  This  paper  is  concluded  with  a 
discussion  of  the  results  and  future  work. 

Design  Problem  and  Analysis 

In  this  section,  the  lumped  element  model  of  a 
piezoelectric-driven  synthetic  jet  is  reviewed  and  the 
basic  dynamic  behavior  discussed. 

Lumped  Element  Model 

At  low  frequencies,  where  the  characteristic 
length  scales  of  the  governing  physical  phenomena 
are  much  larger  than  the  largest  geometric  dimension, 
the  governing  partial  differential  equations  of  the 
dynamic  system  can  be  “lumped”  into  a  set  of 
coupled  ordinary  differential  equations.10  The 
resulting  lumped-parameter  system  can  then  be 
represented  as  an  equivalent  electrical  circuit 
possessing  idealized  discrete  circuit  elements  and 
conjugate  power  variables  for  the  equivalent  voltage 
and  current.11  This  approach  provides  a  simple 
method  to  estimate  the  non-linear  dynamic  response 
of  a  synthetic  jet  actuator  for  design  and  optimization 
purposes. 

A  cross-sectional  schematic  and  corresponding 
equivalent  circuit  representation  for  a  typical 
piezoelectric-driven  synthetic  jet  are  shown  in  Figure 
1.  The  details  of  the  circuit  parameter  estimation 
techniques,  assumptions,  and  limitations  are  given  by 
Gallas  et  al.8  The  structure  of  the  equivalent  circuit  is 
explained  as  follows.  A  harmonic  voltage  Vac  is 
applied  across  the  piezoceramic  to  create  an  effective 
acoustic  pressure  that  drives  the  diaphragm  into 
motion.  This  represents  a  conversion  from  electrical 
energy  to  acoustic  energy  and  is  represented  by  an 
ideal  transformer  possessing  a  turns  ratio  (j)a .  The 
motion  of  the  diaphragm  (i.e.,  volume  velocity,  Q ) 
can  store  potential  energy  via  compressibility  effects 
in  the  cavity  ( Qc )  and/or  can  store  kinetic  energy  via 
oscillatory  flow  through  the  orifice  ( Qout ). 
Physically,  this  is  represented  as  a  volume  velocity 
divider,  Q  =  QC+  Qoul . 
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There  are  several  simplifying  assumptions 
employed  in  the  current  model.  First,  the  synthetic 
jet  is  assumed  to  exhaust  into  a  semi-infinite 
quiescent  air  medium.  In  practice,  these  devices 
interact  with  a  boundary  layer  that  greatly  alters  the 
jet-exit  velocity  profile  and  thus  the  total  orifice 
impedance.12  The  corresponding  differences  in  the 
frequency  response  functions  of  a  synthetic  jet 
actuator  exhausting  into  a  quiescent  medium  versus 
one  exhausting  into  a  boundary  layer  means  that  the 
bench-top  calibration  of  these  devices  is  insufficient 
to  accurately  estimate  the  volume  velocity  or 
momentum  flux  exhausted  into  a  cross  flow  for  a 
given  excitation  voltage.  Second,  compressibility 
effects  in  the  orifice,  but  not  in  the  cavity,  are 
neglected.  The  incorporation  of  these  effects  into  the 
lumped  element  model  is  an  ongoing  research  area 
both  in  the  flow  control  actuator  community  as  well 
as  the  engine  nacelle  acoustic  liner  community.13 

Equivalent  Circuit  Model 

Before  conducting  a  formal  optimization,  the 
basic  dynamic  features  of  the  system  are  reviewed. 
The  frequency  response  function  of  the  volume  flow 
rate  through  the  orifice  per  applied  voltage  for  the 
equivalent  circuit  shown  in  Figure  1  is  a  four-pole, 
single-zero  dynamic  system,8 


Qom  (5)  _ _ _ 

Vac  (s)  a4s4  +  a3s3  +  a2s2  +  axs  +  1 


{1} 


where 

a\  ~  CaD  (Ra0  +  RaN  +  R(iD  ) 

+^aC  {RaO  +  RaN  )  * 

a2  ~  CaD  aRad  ^ aN  ^ aD  ) 

+^aC  {^aRad  +  ^ aN  )  +  ^aC^aD^aD  (^aO  +  RaN  )  >  {2} 


a3 

a4 


:  C  C 

^aC^aD 


MaD  ( RaO  +  RaN  ) 

+  {MaRad  +  M aN  )  RaL 

^aC^aD^ aD  aRad  +  ^ aN  )  * 


,  and 


To  first  order,  the  coefficients  in  Eq.  {2}  are 
constants  determined  via  simple  algebraic 
expressions  as  a  function  of  geometry  and  material 
properties.  This  model  includes  one  empirical 
constant,  RaD ,  that  represents  the  structural  damping 
of  the  piezoelectric  composite.  In  general,  some  of 
the  coefficients  exhibit  frequency  and  amplitude 
dependence  (i.e.,  due  to  nonlinear  effects). 
Specifically,  Ra0  is  a  non-linear  orifice  resistance 

that  is  proportional  to  the  volume  velocity,  Qout  (s) . 

For  a  dc  voltage  ( s  =  0 ),  the  volume  velocity  is 
zero.  At  low  frequencies  ( ^  — >  0  ),  the  volume 


velocity  is  proportional  to  jcodaVac .  At  high 
frequencies  ( s  — »  °°  ), 


Kc  CaCCaDM  aD  (MaN  +MaRad)s3 

and  the  output  attenuates  at  60  dB/decade. 

The  four-pole  system  in  Eq.  {1}  possesses  two 
resonant  frequencies,  j\  and  f2  >  fx ,  that  are  related 
to  the  short-circuit  piezoelectric  diaphragm  natural 
frequency  f D  , 


/o=  — 


27t\lMaDCaD 


and  the  Helmholtz  resonator  frequency  fH  , 


f„=— 


2n\(MaN+MaRad)CaC  ’ 


{4} 


{5} 


by  the  equality 


fj2  =  fofn  ■  {6} 

From  Eq.  {6},  it  is  obvious  that  the  cavity 
volume  and  the  orifice  dimensions,  as  well  as  the 
piezoelectric-diaphragm  characteristics  determine  the 
dynamic  response  of  the  synthetic  jet. 

Model  Verification 

Gallas  et  al.8  have  experimentally  validated  the 
lumped  element  model  for  two  different  prototypical 
synthetic  jet  actuators  using  phase-locked  Laser- 
Doppler  Velocimetry.  For  illustration  purposes, 
these  results  are  briefly  reviewed.  The  dimensions 
and  properties  of  the  devices  are  given  in  Table  1  and 
Table  2.  The  amplitude  of  the  piezoelectric 
excitation  voltage  was  25  V  in  all  cases.  The 
comparison  between  the  full  nonlinear  model 
prediction  and  the  experiments  are  shown  in  Figure  2 
and  Figure  3.  In  both  figures,  the  centerline  velocity 
magnitude  is  plotted  as  a  function  of  frequency.  For 
the  lumped  element  model  predictions,  the  centerline 
velocity  was  estimated  by  modeling  the  flow  in  the 
orifice  as  flow  in  a  circular  duct  driven  by  an 
oscillatory  pressure  gradient.  The  first  case  (I,  Figure 
2)  clearly  illustrates  the  two  resonant  frequencies  of 
the  coupled  oscillator.  The  second  case  (II,  Figure  3) 
corresponds  to  a  system  possessing  a  single  dominant 
peak.  The  lower  peak  at  /,  ~  315  Hz  is  heavily 
damped  in  this  case  due  to  the  frequency  dependent 
nonlinear  orifice  resistance  term  Ra0  .  In  both  cases, 
there  is  sufficient  agreement  between  prediction  and 
experiment  to  justify  the  employment  of  FEM  as  a 
design  and  optimization  tool.  It  is  important  to  note 
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that  the  underlying  assumptions  used  to  derive  each 
lumped  element  limit  the  applicable  frequency  range 
of  the  corresponding  element  from  dc  to  some  upper 
limiting  frequency  (see,  for  example,  Rossi10). 


Table  1:  Piezoceramic  diaphragm  details. 


Case 

Shim  (Brass) 

I 

II 

Elastic  Modulus  (Pa) 

8.963X1010 

Poisson’s  Ratio 

0.324 

Density  (kg/m3) 

8700 

Thickness  (mm) 

0.20 

0.10 

Diameter  (mm) 

23.5 

37 

Piezoceramic  (PZT-5A)  | 

Elastic  Modulus  (Pa) 

6.3xl010 

Poisson’s  Ratio 

0.33 

Density  (kg/m3) 

7700 

Thickness  (mm) 

0.11 

0.12 

Diameter  (mm) 

20.0 

25.0 

Rel.  Dielectric  Constant 

1750 

d3l  (m/V) 

-1.75X1010 

CEF  (F) 

4.42  xlO'8 

6.33xl0‘8 

Table  2:  Synthetic  jet  details. 


|  Case 

Cavity: 

I 

II 

Volume  V0  (m3) 

2.50xl0‘6 

5.00xl0‘6 

Orifice: 

Radius  a0  (mm) 

0.825 

0.42 

Length  L  (mm) 

1.65 

0.84 

Optimization  of  a  Synthetic  Jet 

The  lumped  element  model  presented  in  the 
previous  section  is  a  powerful  design  tool  that 
enables  the  multi-energy  domain  dynamic  modeling 
of  a  synthetic  jet  actuator.  In  this  section,  the  model 
is  used  as  a  vehicle  to  enable  optimization.  For 
simplicity  and  physical  insight,  the  optimization 
problem  is  decoupled  into  two  parts.  First,  the 
piezoelectric  diaphragm  is  held  constant,  and  the 
cavity  volume  and  orifice  dimensions  are  varied. 
Two  different  cost-functions  and  sets  of  constraints 
are  explored  for  this  problem.  Next,  the  piezoelectric 
diaphragm  is  optimized  and  various  concepts  for  the 
driver  configuration  are  presented  along  with  their 
associated  results.  In  all  cases,  the  optimization 
problem  was  solved  via  MATLAB’s  optimization 
toolbox. 

Orifice/Cavity  Optimization 

For  a  given  piezoelectric  composite  diaphragm, 
the  device  behavior  is  governed  by  the  cavity  volume 
V0 ,  the  orifice  radius  aQ ,  and  orifice  length  L  .  For 
example,  Figure  4-Figure  6  illustrate  the  effect  of 


varying  V0 ,  a0 ,  and  L  on  the  centerline  velocity 
output  of  the  “nominal”  synthetic  jet  actuator  shown 
in  Figure  2  (Case  I).  Increasing  the  cavity  volume 
results  in  a  larger  acoustic  compliance,  CaC  ,  and  thus 

lowers  first  resonance,  /j  oc  \j  ,  and  also 

decreases  the  overall  broadband  amplitude,  while 
slightly  reducing  the  second  resonant  frequency 
relative  to  the  nominal  case  (Figure  4).  Increasing 
the  orifice  radius  decreases  total  acoustic  mass, 
MaN  +  MaRad  ,  resulting  in  a  higher  first  resonance 
and  a  correspondingly  higher  second  resonance 
(Figure  5).  Conversely,  increasing  the  orifice  length 
increases  the  total  acoustic  mass,  MaN  +  M aRad  ,  and 
thus  reduces  the  first  resonance,  much  in  the  same 
manner  as  varying  the  cavity  volume  (Figure  6).  It  is 
clear  that  varying  these  three  geometric  variables  can 
significantly  affect  the  frequency  response  function. 

Different  flow  control  applications  will  require 
different  frequency  response  specifications  in  terms 
of  bandwidth  and  the  output  physical  quantity  of 
interest.  For  example,  certain  applications  may 
require  significant  actuation  authority  in  terms  of 
momentum  flux  over  a  narrow  frequency  range, 
while  another  may  require  a  flat,  broadband  response 
possessing  a  prescribed  minimum  centerline  velocity. 
This  emphasizes  the  importance  of  clearly  specifying 
a  design  objective. 

In  addition  to  choosing  a  design  objective,  the 
optimization  problem  requires  the  specification  of 
various  constraints.  These  constraints  can  be  broadly 
classified  into  three  categories:  performance 
constraints,  design  variable  constraints,  and  model 
constraints.  Performance  constraints  may  consist  of 
specifying  the  desired  locations  of  the  resonant 
frequencies,  the  flatness  tolerance  of  the  gain  factor 
of  the  frequency  response  function,  etc.  Design 
variable  constraints  may  be  needed  due  to  physical 
limitations  based  on  manufacturability,  packaging 
considerations,  material  failure,  etc.  As  stated  at  the 
end  of  the  previous  section,  there  are  limits  to  the 
applicability  of  each  lumped  element  as  a  function  of 
frequency.  Therefore,  after  an  optimal  solution  is 
achieved,  the  validity  of  the  assumptions  in  the  model 
must  be  evaluated,  or  the  optimization  problem  must 
be  constrained  by  the  limitations  of  the  model.  The 
latter  approach  was  chosen  for  this  study.  For 
example,  the  lumped  elements  for  the  piezoelectric 
composite  are  limited  to  frequencies  below  the 
second  natural  frequency  of  the  composite 
diaphragm.  The  size  of  the  cavity  is  limited  such  that 
the  impedance  of  the  cavity  can  be  approximated  by  a 
compliance.8  Another  limitation  appears  for  the 
orifice  aspect  ratio,  L/a0  .  Based  on  experimental 
results  for  several  cases,  Gallas9  found  that 
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reasonable  agreement  was  achieved  between  the 
lumped  element  model  and  measured  dynamic 
response  when  the  orifice  aspect  ratio  approximately 
exceeded  unity. 

For  the  present  optimization  study,  the  goal  is  to 
improve  the  performance  of  the  nominal  designs 
presented  in  Figure  2  (Case  I)  and  Figure  3  (Case  II). 
The  objective  function,  constraints,  and  design 
variables  are  summarized  in  Table  3.  Specifically, 
Case  I  is  optimized  using  one  cost  function,  while 
Case  II  is  optimized  using  two  different  cost 
functions. 

The  first  cost  function  employed  maximizes  the 
integrated  centerline  velocity  over  the  entire 

frequency  range,  [/hm u0(f)df  where  the  upper  limit 

of  integration  is  3000  Hz  and  1500  Hz  for  Case  I  and 
II,  respectively.  The  motivation  for  such  an  objective 
function  is  to  increase  the  broadband  response  of  the 
actuator.  The  constraints  on  the  orifice  radius  are 
motivated  by  device  manufacturability  and  flow 
perturbation  concerns,  while  the  minimum  orifice 
length  constraint  is  driven  by  the  requirement  that  the 
orifice  plate  be  rigid.  The  volume  range  is  dictated 
by  size  limitations.  In  addition,  constraints  are 
imposed  on  the  orifice  aspect  ratio  and  by  fixing  the 
piezoelectric  drive  voltage.  Figure  7  shows  the 
resulting  optimized  frequency  response  function 
compared  to  the  nominal  response  for  Case  I.  The 
frequency  response  function  increased  over  the  entire 
frequency  range  by  decreasing  V0 ,  a0,  and  L .  In 
addition,  the  first  resonant  peak  is  heavily  damped  in 
the  optimized  response. 

The  results  using  a  similar  approach  for  Case  II 
are  shown  in  Figure  8.  Again,  the  frequency 
response  function  increased  over  nearly  the  entire 
frequency  range  by  decreasing  V0 ,  a0 ,  and  L .  In 
addition,  the  flat  portion  of  the  response  function  is 
increased  and  the  second  resonant  peak  is  moved  to  a 
higher  frequency.  This  actuator  design  is  useful  for 
applications  that  require  a  flat  broadband  response. 

The  second  cost  function  employed  maximizes 
the  centerline  velocity  at  the  second  resonant 
frequency  of  the  system,  u0(f2)  .  From  a  practical 
standpoint,  this  optimization  is  useful  for  applications 
requiring  high  actuation  authority  over  a  narrow 
frequency  range.  The  constraints  on  the  orifice  and 
cavity  geometry,  as  well  as  the  aspect  ratio  and  drive 
voltage  are  similar  to  those  outlined  above.  A  new 
equality  constraint,  however,  is  placed  on  the  location 
of  the  second  resonant  frequency.  The  results  of  this 
optimization  for  Case  II  are  shown  in  Figure  9.  The 
second  resonant  peak  is  increased  as  desired  (by  ~ 
35%).  In  addition,  the  broadband  response, 


especially  at  the  lower  frequencies,  is  increased  with 
respect  to  the  nominal  case. 

While  the  optimization  studies  in  this  section 
indicate  the  promise  of  improved  performance, 
additional  gains  can  likely  be  achieved  by  optimizing 
the  piezoelectric  composite  diaphragm.  This  is 
described  in  the  next  section. 


Table  3:  Summary  of  optimization  problem  for  a 
_  synthetic  jet  actuator.  _ 


Objective 

Constraints 

Variables 

J0/l,m  uo  (/)# 

•  Upper/Lower  bounds  on 
variables: 

Case  I: 

a0 

L 

Vo 

f  0.125  <  ao  <2 
\  0.5  <  L  <  20 
[ll52<  V0  <47,640 

Case  II: 

f  0.125  <  aQ  <2 

0.38  <  L  <  20 
[780  <Vo  <21,270 

•  Orifice  aspect  ratio 
L/a^  >  1 

•  Fixed  input  voltage 
V  =25V 

ac 

Uo(fo) 

•  Upper/Lower  bounds  on 
variables: 

Case  II: 

a0 

L 

Vo 

f  0.125  <ao  <2 

0.38  <  L  <  20 
[780  <y  <21,270 

•  Natural  frequency 
f2  =  820  Hz 

•  Orifice  aspect  ratio 
L/a^  >  1 

•  Fixed  input  voltage 
y  =25V 

ac 

Piezoelectric  Composite  Optimization 

For  a  given  cavity/orifice  configuration,  the 
response  of  the  synthetic  jet  is  directly  proportional 
to  the  effective  acoustic  piezoelectric  coefficient,  da  , 
as  shown  in  Eq.  { 1 } .  The  response  is  also  governed 
by  the  natural  frequency  of  the  diaphragm,  fD , 
which  is  a  function  of  the  acoustic  mass,  M aD  ,  and 
short-circuit  acoustic  compliance,  CaD ,  of  the 
piezoelectric  diaphragm.  These  lumped  elements  are 
all  dictated  by  the  geometry  and  material  properties 
of  the  diaphragm. 

Before  conducting  a  formal  optimization  of  the 
piezoelectric  composite  diaphragm,  the  basic 
modeling  concepts  are  reviewed.  Top- view  and 
cross-sectional  schematics  of  two  axisymmetric 
piezoelectric  composite  plate  configurations 
examined  in  this  study  are  shown  in  Figure  10-Figure 
13.  The  first  type  (Figure  10  and  Figure  11)  consists 
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of  an  inner  circular  disc,  possessing  a  piezoceramic 
material  of  thickness  t  and  radius  Rx  bonded  to  a 

shim  material  of  radius  R2  and  thickness  ts .  The 
second  type  (Figure  12  and  Figure  13)  consists  of  an 
annular  ring  possessing  a  piezoceramic  material  of 
radial  extent  R2  -  Rx  bonded  to  a  shim  material  of 
radius  R2 .  In  general,  the  actuator  is  subjected  to  an 
applied  voltage  Vac  across  the  piezoceramic 
thickness.  This  loading  creates  both  transverse  w(r) 
and  radial  u(r)  displacements. 

Ideally,  a  piezoelectric  actuator  is  a  linear, 
conservative,  reciprocal  transducer.10  The 
piezoelectric  composite  deforms  in  response  to  both 
an  applied  ac  voltage  and  a  differential  pressure.  The 
lumped  piezoelectric  coupling  equations  for  the 
transduction  model  are14 


where  AV  is  the  volume  displaced  by  the  plate  due 
to  the  application  of  a  differential  pressure  P  and/or 
voltage  Vac ,  q  is  the  charge  stored  on  the 
piezoelectric  electrodes,  CEF  is  the  electrical  free 
capacitance  of  the  piezoelectric  material,  CaD  is  the 
short-circuit  acoustic  compliance  of  the  plate,  and  da 
is  the  effective  acoustic  piezoelectric  coefficient. 

The  acoustic  mass  MaD  is  determined  by 
equating  the  lumped  kinetic  energy  of  the  vibrating 
diaphragm  expressed  in  acoustic  conjugate  power 
variables  to  the  total  kinetic  energy.  The 
determination  of  the  lumped  element  parameters  for 
this  two-port  model  requires  the  solution  of  the 
transverse  static  deflection  field  as  a  function  of 
pressure  and  voltage  loading.  For  this  study,  linear 
laminated  plate  theory  is  used  to  solve  for  the 
transverse  static  deflection,  and  an  optimization 
scheme  is  then  implemented  using  this  two-port 
model.  The  details  of  the  composite  plate  model  are 
presented  in  Prasad  et  al.14 

The  material  properties  are  Es ,  vs ,  Ep  ,  and  Vp  , 
where  the  subscripts  4 s’  and  ‘p’  denote  the  shim  and 
piezoelectric,  respectively,  while  E ,  v ,  and  d3l  are 
the  elastic  modulus,  Poisson’s  ratio,  and  piezoelectric 
constant,  respectively.  In  this  study,  the  material 
properties  are  constant  and  are  listed  in  Table  1. 

The  optimization  objective  function  maximizes 
the  achievable  displaced  volume  per  applied  voltage, 
da  =  AV| P=0/Vac  max  ,  compared  to  the  nominal  values 

for  Cases  I  and  II.  The  motivation  behind  this 
objective  function  is  that  the  strength  of  a 


piezoelectric-driven  synthetic  jet  is  directly 
proportional  to  this  quantity.  Note  that  Vacmax  is  the 

maximum  voltage  that  can  be  applied  to  the 
piezoelectric  without  depolarization  and  is  estimated 
as 

Kc.max  =?px£'/max  =tpx(30V/mil),  {8} 

where  Efmax  is  the  coercive  electric  field  applied 
through  the  thickness  t  of  the  piezoceramic. 

Next,  constraints  are  defined.  Reasonable  lower 
and  upper  bounds  are  defined  for  the  geometric 
variables  ts ,  tp  and  Rx ,  while  R2  is  fixed  in  this 

study.  The  resonant  frequency  of  the  composite 
diaphragm  plays  an  important  role  in  the  design 
process.  Without  a  constraint,  the  optimization 
results  in  a  very  compliant  diaphragm  that  can 
achieve  large  volume  displacements  but  has  a  low 
natural  frequency.  Since  bandwidth  is  often  an 
important  quantity,  a  lower  bound  is  placed  on  the 
resonant  frequency  of  the  diaphragm.  The  goal  here 
is  to  improve  the  performance  of  the  diaphragm 
compared  to  the  existing  designs  in  Table  1.  Thus, 
the  fundamental  frequencies  of  the  nominal  designs 
used  in  Cases  I  and  II  are  chosen  as  the  lower  bound. 

Optimizations  of  both  the  inner-disc  and  outer¬ 
ring  configurations  were  thus  performed  for  Cases  I 
and  II  using  the  MATLAB  optimization  toolbox. 
The  results  are  summarized  in  Table  4  for  the  inner- 
disc  configuration.  The  optimum  design  slightly 
improves  (by  ~6%)  the  volume  displacement 
achieved  for  the  nominal  Case  I,  while  a  52% 
improvement  is  achieved  compare  to  the  nominal 
Case  II.  It  should  be  noted  that  the  lower  frequency 
bound  on  f D  is  always  active  in  the  optimized 
design,  indicating  that  a  direct  comparison  of  da 
between  the  original  and  optimized  designs  is 
justified. 


Table  4:  Optimization  results  of  piezoelectric 
diaphragms  with  an  inner-disc  configuration. 


I 

H 

Original 

Optimum 

Original 

Optimum 

R2  (mm) 

11.5 

11.5 

18.5 

18.5 

R}  (mm) 

10.0 

10.5 

12.5 

16.9 

ts  (mm) 

0.15 

0.10 

0.10 

0.08 

tp  (mm) 

0.08 

0.16 

0.11 

0.12 

fD  (Hz) 

2114 

2114 

632 

632 

da  (xio-10) 
(m3/V) 

0.553 

0.586 

4.32 

6.56 

The  optimization  results  are  summarized  in 
Table  5  for  the  outer-ring  configuration.  The 
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performance  improvement  is  3%  and  46%  compared 
to  the  nominal  Case  I  and  II  designs,  respectively. 
The  performance  improvement  is  less  than  the  inner- 
disc  configuration.  However,  the  outer-ring 
configuration  is  less  apt  to  lead  wire  failures  than  the 
inner-disc  case  due  to  the  reduced  motion  near  the 
clamp  (see  Figure  13). 

Table  5:  Optimization  results  of  piezoelectric 

diaphragms  with  an  outer-ring  configuration. 


I 

II 

R2  (mm) 

11.5 

18.5 

R{  (mm) 

7.0 

10.9 

ts  (mm) 

0.13 

0.095 

tp  (mm) 

0.056 

0.048 

f D  (Hz) 

2114 

632 

da  (xlO'10) 
(m3/V) 

0.568 

6.30 

Other  configurations  are  possible  and  will  be 
studied  in  future  work.  For  example,  oppositely- 
polarized  inner-  and  outer-rings  can  be  combined  in 
both  unimorph  and  bimorph  designs. 

Conclusions  and  Future  Work 

The  optimization  of  a  piezoelectric-driven 
synthetic  jet  actuator  based  on  LEM  has  been  carried 
out.  It  has  been  shown  that  LEM  is  a  viable  tool  to 
optimally  design  these  devices  for  candidate 
applications.  To  simplify  the  problem,  the  current 
study  has  split  the  optimization  problem  into  two 
parts  by  separately  optimizing  the  (1)  cavity  volume 
and  orifice  dimensions  and  (2)  the  piezoelectric 
diaphragm. 

A  constrained  optimization  of  the  cavity  volume 
and  orifice  dimensions  of  two  baseline  synthetic  jets, 
each  with  a  given  piezoelectric  diaphragm,  has  been 
conducted  using  two  different  objective  functions. 
One  function  seeks  to  improve  the  centerline  output 
velocity,  which  is  a  representative  measure  of  the 
device  output,  over  a  broad  frequency  range  of 
interest.  The  other  objective  function  focuses  on 
maximizing  the  centerline  velocity  at  a  prescribed 
resonant  frequency  of  the  device.  Significant 
improvements  have  been  achieved  for  all  cases. 

In  addition,  two  candidate  piezoelectric 
diaphragm  configurations  have  been  optimized.  One 
uses  the  standard  inner-disc  piezoceramic  patch 
bonded  to  a  metal  shim,  while  the  other  uses  an  outer 
piezoceramic  ring.  In  each  case,  the  objective  is  to 
maximize  the  achievable  volume  displacement  of  the 
diaphragm,  corresponding  to  the  coercive  electric 
field  strength  of  the  piezoceramic.  The  natural 
frequency  of  the  piezoelectric  diaphragm  was 
constrained  to  be  greater  than  or  equal  to  that  of  the 


baseline  designs.  Both  methods  yield  modest  (~5%) 
improvements  for  one  diaphragm  and  significant 
improvements  for  the  other  diaphragm  (>50%).  The 
inner-disc  configuration  yields  slightly  better 
performance,  but  is  likely  to  be  more  susceptible  to 
lead- wire  breakage. 

In  future  work,  experiments  will  be  conducted  to 
confirm  the  optimization  results  presented  here. 
Other  piezoceramic  diaphragm  configurations  will  be 
studied.  In  addition,  a  coupled  constrained 
optimization  will  be  carried  out  that  seeks  to  optimize 
the  device  performance  by  simultaneously  varying 
the  piezoelectric  diaphragm  dimensions,  cavity 
volume,  and  orifice  dimensions. 

Also,  another  optimization  problem  of  a 
synthetic  jet  will  be  investigated.  The  current  study 
is  limited  to  improving  an  existing  baseline  design. 
Perhaps  the  more  interesting  case  is  the  optimal 
design  synthesis  problem.  In  this  problem,  the 
designer  seeks  to  achieve  a  desired  frequency 
response  function.  Due  to  the  nonlinear  nature  of  the 
system,  the  design  objective  can  be  approximated  by 
a  linear  transfer  function  that  is  valid  at  a  particular 
driving  voltage.  A  key  challenge  here  is  that  the  end 
user  must  be  able  to  translate  desirable  actuator 
characteristics  into  quantitative  design  goals. 
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domain  domain 

Figure  1:  Schematic  representation  and  equivalent 
circuit  model  of  a  piezoelectric-driven  synthetic  jet. 


Figure  2:  Comparison  between  the  lumped  element 
model  and  experiment  for  Case  I.8 
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Figure  3:  Comparison  between  the  lumped  element 
model  and  experiment  for  Case  II.8 


Figure  4:  Lumped  element  model  prediction  for 
varying  cavity  volume  ( Case  I). 


Figure  5:  Lumped  element  model  prediction  for 
varying  orifice  radius  ( Case  I). 


Figure  6:  Lumped  element  model  prediction  for 
varying  orifice  length  ( Case  I). 


Figure  7:  Optimization  of  Case  I,  maximizing  the 
overall  output. 


Figure  8:  Optimization  of  Case  II,  maximizing  the 
overall  output. 
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Figure  9:  Optimization  of  Case  II,  maximizing  the 
second  peak. 


Figure  10:  Top  view  of  the  inner-disc  piezoelectric 
composite  diaphragm  configuration. 


Figure  12:  Top  view  of  the  outer-ring  composite 
piezoelectric  diaphragm  configuration. 


Figure  13:  Cross-section  of  the  outer-ring 
piezoelectric  composite  diaphragm  configuration. 
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Figure  11:  Cross-section  of  the  inner-disc 
piezoelectric  composite  diaphragm  configuration. 
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